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PROGRAM SUMMARY 



• Title of program:. SANC 

• Catalogue identifier. 

• Program obtainable from:. Internet sites at DLNP. JINR, Dubna, Russia, 
http://sanc.jinr.ru/ (159.93.74.10) and at CERN, http://pcphsanc.cern.ch (137.138.180.42) 

• Computers for which the program is designed and others on which it has been tested: 
Designed for: platforms on which Java and FORMS are available 

Tested on: Intel-based PC's 

• Operating systems: Linux, Windows 

• Programming languages used: Java, FORMS, PERL, FORTRAN 

• Memory required to execute with typical data: 10 Mb 

• No. of hits in a word: S2 

• No. of processors used: 1 on SANC server, 1 on SANC client 

• Distribution form,a,t\ gzipped tar archive 

• Keywords: Feynman diagrams. Perturbation theory. Quantum field theory. Standard Model, Elec- 
troweak interactions, QCD, QED, One-loop calculations, Monte Carlo generators. 

• Nature of physical problem: Automatic calculation of pseudo- and realistic observables for various 
processes and decays in the Standard Model of Electroweak interactions, QCD and QED at one- 
loop precision level. Form factors and helicity amplitudes free of UV divergences are produced. For 
exclusion of IR singularities the soft photon emission is included. 

• Method of Solution: Numerical computation of analytical formulae of form factors and helicity 
amplitudes. For simulation of two fermion radiative decays of Standard Model bosons {W"^, Z) and 
the Higgs boson a Monte Carlo technique is used. 

• Restrictions on the complexity: In the current version of SANC there are 3 and 4 particle processes 

and decays available at one-loop precision level. 

• Typical Running time: The running time depends on the selected process. For instance, the symbolic 
calculation of form factors (with precomputed building blocks) of Bhabha scattering in the Standard 
Model takes about 15 sec, helicity amplitudes — about 30 sec, and bremsstrahlung — 10 sec. The 
numerical computation of cross-section for this process takes about 5 sec (CPU SGHz IP4, RAM 
512Mb, L2 1024 KB). 
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1 Introduction 



Project motivation 

The main goal is tlie creation of a computer system for semi-automatic calculations of realistic and 
pseudo-observables for various processes of elementary particle interactions "from the SM Lagrangian to 
event distributions" at the one-loop precision level for the present and future colliders - TEVATRON 
(Runs II and III), LHC, electron Linear Colliders (ISCLC, CLIC), muon factories and others. 

Furthermore, the SANC system, even at the level which is reached already, may be used for educational 
purposes by students specializing in high energy physics. With its help, it is easy to follow all steps of 
calculations at the one-loop precision level for W,Z,H^ ffil), H 77, Z7, ZZ, WW, t bW decays, 
and many other processes. Moreover, all the calculations arc realized in the spirit of the book [1] which 
makes the SANC system particularly appealing for pedagogical purposes. 

Historical overview 

The SANC project has been started in early 2001. During the first phase of the project (2001-2003), 
the SANC group demonstrated the workability of the computer system which is being developed [2] . The 
version 0.01, from 03/28/2001, was already able to compute one-loop Feynman diagrams for all SM 
1 — > 2 decays and 2/ — > 2/ processes (in i?| and unitary gauges, including QCD, accessing thereby 
all one- loop diagrams needed for the processes considered by the Dubna group in the past [3]- [6] in 
connection with theoretical support of experiments at CERN and DESY. The FORM codes (at present 
FORMS [7] is being used) computing the one-loop ultraviolet finite scalar form factors of the amplitudes 
of the decays Z{H, W) ff were unified and put into a special program environment, written in JAVA. 
This version was used for a revision of Atomic Parity Violation [8] , and for a calculation of the one- loop 
electroweak radiative corrections for the process e+e~ — > // [9] and neutrino DIS [10]. 

In the second phase of the project (2004 2006), we extend automatic calculations of such a kind to a 
large number of HEP processes, with emphasis on LHC physics. 

Present status 

The present level of the system is realized in the version vl .00. This version has a fresh new layout 
and is more user friendly than earlier versions. 

New in version 1 . 00 are Compton scattering and several other ffbb processes. By our philosophy we 
treat them as building blocks for future calculations of 5 ^ processes (fully massive case). 

For the last year we substantially enhanced our computer system compared to the status presented 
in the years 2002-2003 at large-scale international conferences, such as ACAT2002 [11]-[14], ICHEP2002 
[15], RADCOR2002 [16] and Workshops at Saint-Malo [17], CERN [18], Montpelier [19] and Paris [20]. 

How to get started and use SANC 

To learn more about available SANC servers look at our home pages at Dubna http://sanc.jinr.ru or 
CERN http://pcphsanc.cem.ch. 

SANC may be accessed via the so called SANC client a software free to download. The user will 
always get the latest (updates) versions from either of the two above addresses. 

Levels of the calculations 

SANC is subdivided into three logical levels, each with a specific purpose. 
• Level 1, Analytic 

The analytical application includes enhanced tools of FORM procedures. SANC has three types of 
procedures: specific, intrinsic and special. The specific procedures are used by FORM source codes, typically 
only once; they are always visible (see Section 3) and can be modified by the user. The action of intrinsic 
procedures is uniquely specified by their arguments, therefore they may be used by many FORM codes. 
Their bodies are not accessible to the users. Finally, special procedures are used only a few number of 
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times to perform some special action in a given FORM code. Normally, they have very simple arguments 
like field indices. 

The Figures 2 3 show fully open menus for "Precomputation" and for available "Processes" in the 
QED part, Fig f8, and "Processes" in the EW part. Fig. 19. In this article we explain in detail the 
process of "Precomputation" . 

Entering your chosen process, you are in an active session and receive the analytical result for scalar 
Form Factors (FF), Helicity Amplitudes (HA) and the accompanying Bremsstrahlung contributions (BR). 

For the calculation of the HAs we use techniques of Vega-Wudka [21] and Kleiss-Stirling [22]. 

As a main example of the description of the calculation of FF and HA for processes /i/i ^ // we 
mention Ref. [9] . All calculations of FF^HA^BR on the SANC tree are realized in the same job stream. 

• Level 2, Numerical 

The analytic results are transferred to the second level where they are analyzed by a software package 
s2n.f {sym,bols-to-num,bers), written in PERL. The s2n.f package automatically generates FORTRAN codes 
for subsequent numerical computations of decay rates and process cross sections. The calculational flow 
inside levels 1 and 2 and the exchange of data between them is fully automatized and is governed by 
selecting corresponding items in menus. 

The s2ii part of the SANC system is completely fixed for all available decays (besides {Z W^W~); for 
NC 4/ — > processes we have s2n for FF— >HA— >SoftBR parts, and for CC processes fif[ — > we have 
s2n for FF^HA^ SoftBR+HardBR parts (here /i with or without primes denotes massless fcrmions, 
e.g. of the 1st generation). We have performed many high-precision comparisons of the numerical results 
derived with the aid of s2n.f with the results of the alternative systems FeynArts [23] and GRACE [24] 
and the code topf it [25]. 

• Level 3, MC generators 

In version 1 . 00, MC generators are available only for B ff decays together with relevant graphic in- 
terface. The results can be presented in a variety of histograms. The user may "play with the parameters" 
of histograms in the window. 

A first Monte Carlo generator for decays Z{H, W) — » //7 is created in close contact with members 
of the KK collaboration, see papers [26, 27]. The MC generator is also accessible via menus, therefore, 
for the case of decays we are able to demonstrate how the full cJiain of ealculations works out within our 
integrated system. The MC event generators are supposed to be usable also in a "stand alone" mode 
ready to be incorporated into the software of experiments. 

This paper is organized as follows: 

In Section 2 we describe amplitudes for all available in version 1.00 3- and 4-leg processes. 

Section 3, the main section of this paper, is fully devoted to the concept of precomputation; a com- 
prehensive description of the SANC precomputation tree and its modules is given. We assume that while 
reading this part the reader will be looking inside the corresponding modules. For one example, namely 
precomputation of photonic vacuum polarization, we demonstrate the whole process of calculations by 
presenting intermediate results after each procedure is called. 

In Section 4 we describe a part of SANC procedures, mostly those which are used by precomputation 
modules. 

In Section 5 we briefly describe the SANC trees of processes implemented for the time being. 

Although this paper is mostly devoted to the SANC precomputation, in Section 6 we give a short User 
Guide for version 1 . 00. A more detailed description of the SANC processes trees and computer aspects 
of SANC system will be given elsewhere. 
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2 Amplitude Basis, Scalar Form Factors, Helicity Amplitudes 



2.1 Introduction 

In this section we present a collection of formulae for the amplitudes of basic SM 1^2 decays and 2^2 
processes available in SANG v. 1.00. The covariant one-loop amplitude (CA) corresponds to a result of 
the straightforward standard calculation by means of SANG programs and procedures of all diagrams 
contributing to a given process at the tree (Born) and one-loop levels. It is represented in a certain basis, 
made of strings of Dirac matrices and/or external momenta (structures), contracted with polarization 
vectors of vector bosons, if any. We usually omit Dirac spinors. The amplitude also contains kinematical 
factors and coupling constants and is parameterized by a certain number of form factors, which we denote 
by J-', in general with an index labeling the corresponding structure. The number of FFs is equal to the 
number of structures. If there is only one FF, we normally do not label it. For the processes with non 
zero tree-level amplitudes the FFs have the form 

J^=l + k^, (1) 

where "1" is due to the Born level and the term ^ with the factor 

167r2 ' 

is due to the one-loop level. We also use various coupling constants 



k=l^, (2) 



Qf , if^ , af = Vf + ttf , 5f = vf - af , s„ = - , c„ = , etc. (3) 

Given a CA parameterized by a certain number of FFs, SANG computes a set of HAs, denoted by 
WA1A2A3...) where A1A2A3 . . . denote the signs of particle spin projections onto a quantization axis as will 
be explained in the following sections. 

In the representation of massive HAs, the following notation is very useful: 



P±(/,Mi,M2) = •\/|/-(Mi±M2f|. (4) 
2.2 The 3-leg processes B{Q) f{pi) + /(ps) 

In this section we present amplitudes for 1 — > 2 decays involving one vector boson and two fermions. For 
all B ^ ff decays, except Higgs boson decay, the three Aj in HA1A2A3 denote the signs of the boson, 
fermion and antifermion spin projections, respectively. Fermion spins are projected onto their momenta, 
the boson spin is projected onto the fermion momentum. For example, corresponds to the following 

spin configuration: 

^ ^ 

Spin B 



P2 Pi 



Spin / Spin / 
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Its CA is described by one FF only: 



Correspondingly, there is only one independent HA: 

7i++ = W— 
— = TH. — I- = 0, 

here 

P+„ = P+ {Mlmf,mf) = ^Ml-Am). 

The CA of the decay of a heavy vector particle (V) contains two FFs: 

AAff = eQf [i"ffj.TQ + mfDf^To] £^(0) > 

here and below in this section, until stated otherwise, = {p2 — Pijfi- 
The two independent HAs are: 

'Ho++ = Ho — =emfQf {P^^^Y To + 2J^q 
n+±+ = n± — = rto+- = rto-+ = rt_+± = 



( 



P+„ = P+ {Ml,mf,mf)=^M^-Am). 



In this case, the CA is described by three FFs: 



A 



9 
2c„ 



where 76, 7 = 1 ± 75. The three independent HAs look as follows: 



H++ 
H 



+ — 



Hq 



9 

\/2cw 



lf\M^ + P+„)T^ + SfM^T^ 
gnif 

Cw ^ 

= TioH — = ^0 — h = — 1-+ = ^ — 1-- 



P+„ = P+ {Ml rrif, mf) = J - Am} . 



W+ ^u + d (W^ -^ u + d) 

The CA of this decay is described by four FFs: 

9 



The corresponding HAs read: 

^ + = I {^wud + ^wud) ^i- + {-^wud ~ ^wud) ' 

0++ = -TTIkTI- {Pwud'^n+d + P^ud'^u-d)^L + {P„ud'^n+d-P^ud'^u-d)^R 



(14) 



'^^wud^wud [^wud "I" -^wud] "I" ^wud-^wud [-^wud ~ ^wud\ 

9 1 



2^/2 M„ 



{Pwud^u+d - P+^^niu-d) Tr. + (P^„rfm„+d + P+„^m„_d) 



+ P+ (' P+ — P^ P 4- P+ P^ (' P+ -I- P^ P 

^ K'ud wud \ wud wud) ld ^ wud wud \ wud ^ wudj « 

7Y+±+ = 'H±-- = Ho^ — = Ho — I- = H — 1-± = , 



where 



^wud = Pi^w, m„, ±md) = - m^^^ , m„±d = m„ ± ma . 



(15) 
(16) 



t^b + W+ 

The CA has the same structure as in Eq. (14), but here with D^^ = (pi +p2)n- 
The four HAs are: 

9 
2 



n. 



-0 = — 



(^tt6"^*-b - Ptwb^t+b) J^L - (P4+ femt_6 + P^^^mt+b) Tk 

'^P-twbPtwb [Ptwb + Ptwb) -^i^o — Ptwb^twb {^tivb ^ Ptwb) -^^o ' 

{Ptwb'T^t-b + P~^f,mt+b) - {P^^^mt-b - P-^t,mt+b) Tn 



2\/2M„ 



^Ptv/b^tv/b ^Ptwb Ptwb) -^^^ PtwbPtwb ^Ptwb + -^twb) -^"o 

W+_o = W-i — = W — 1-+ = W_+o = W — ± = 0, 



where 



P{M„, rrit ± mfc) = J m^^;, - M^, , mt±b = mt±mb. 



(17) 

(18) 



The amplitudes of this process are similar to the previous one (though not identical). Their exphcit 
expressions can be found in the relevant module of the SANC tree. 
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2.3 The 3-leg processes B{Q) V{pi) + V{p2), V = 'y,Z,W 

In this section we just list CAs and HAs for basic three-boson decays in the SM. Note, the first two 
decays do not proceed at the tree level, this is why their FFs do not start with "1" . 



H 7(pi) + 7(p2) 



Ml 



W++ = H — = kgSy^T, 
n+- = H-+ = 0. 



Ah~,z = ikgsy 
H++ 



Ml 



( M^\ ~ 

n- = kgs^ (i-^j-^, 
Wo± = W — I- = . 



H ^ Z{p^) + Z{p2) 



gM, 

Cw 



Sf^u^Fo + ^-J^^P ] <^iyiPl)<^i^{P2) , 



++ 



Woo 



n. 



gMz 



gM^ 



2 Ml) 

H±Q = W-i — = Wo± = — h = . 

H ^ W{pi) + W{p2) 

•A.HWW 



4M2 



{-gMw) ( S^uJ^D + ^^^i^^p ] eL/(Pi)e^(P2) , 



++ 



n. 



{-gM„)Tn, 



IMl 
2M2 



1 - 



IMl 
4M2 



Woo = {-gM^) 

Ti.+o = H-i — = Wo+ = 'Ho- = Ti. — i- = W_o = . 
Z{Q) ^ Wipi) + W{p2) 

•Azww — g C-w 

+ (f}af,yDpT^°l^ ea(Q)ei.(pi)e^(p2) 
The HAs for this decay are not implemented in SANC v. 1.00. 



D2 



(19) 



(20) 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 
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2.4 The 4-leg NC processes fji (jZ) ff 

Here we present the CAs and HAs for any /i/i// NC process at any channel s, t or u. Here 
stands for vacuum, and by /i we mean a first generation fermion with field index 11,12,13,14, whose mass 
is neglected everywhere except in arguments of logs (mass singularities) and by / we mean any fermion 
with field indices from 11 to 22 (see Section 6 for definition of field indices). For such a case, the Higgs 
and (ffi boson interactions with the fi current are also neglected. 
The covariant one-loop amplitude of the 2/ 2/ process 




Figure 1: /i/i ^ // process. 

is described in terms of six form factors: LL, QL, LQ, QQ, LD and QD, corresponding to six Dirac 
structures {A-y is also described by a QQ structure, it is separated out for convenience; a{s) = a(0)jF^(s)). 
Note that all 4-momenta are incoming and the usual Mandelstam invariants in our metric (i.e. = — m^) 
are: 

{Pl+P2f = -S, {P2+P3f = -t, {P2+P4f = -U. (28) 

The 7 and Z exchange amplitudes are: 

A^^s) = ie2^^.F^(s)7^®7^ , (29) 
s 

X l^l'fflfln (1 + 75) ® 7/. (1 + 75) :Fi.l{s, t, u) +Sfjfj^ 75) :Fql{s, t, u) 

+IffSf-ff, (1 + 75) «) -ffj,J^LQ{s,t,u) + Sf^df-f^ (g) -ff,J^QQ{s,t,u) (30) 

+4f 4^^7(u(l + 75) O {-imfD^)J^LD{s,t,u) + dfjf^^^ (g) {-imfDij,)J^Qo{s,t,u)^. 

Here and below in this and in the next sections, = {p4 — ps)^ and Xz(s) is the Z/7 propagator ratio 

1 s 

Symbol (g 7^ is used in the following short-hand notation: 

1ti®li^ = v (pi) j^u {P2) u ipz) j^v (pa) . (32) 

For more details see Ref. [9]. 
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If the /i mass is neglected, we have six corresponding HAs. They depend on kinematical variables, 
coupling constants and our six scalar form factors: 



++- — 



H-+-+ 



+-++ 



-e2(l + cost?) \ Qf,QfJ^^{s)+Xz{s)6f, 
-e2(l-cost?) (gA0/^^(s)+x4s)% 



n-+++ = ^ sin-ff [Qf^Q (s) 



H+ = 



„2 2m/ 



-e — ^ sin 1? I g/jQ/jP^(s) +Xz(s) 



(3), 



-e2 (1 - COST?) \Qf,Qf:F^{s) + Xz{s) 



H+—+ = -e^{l + cos^)\Qf,QfT^{s)+Xz{s) 
W+++± = W++_± = H — +± = H ± = 0; 



(33) 



helicity indices, for example, H h denote the signs of the fermion spin projections onto their momenta 

Pi,P2,P3,P4, respectively. 
Moreover, 



P"*" = P~^{s,mf,mf) = - 4m J . 



and the scattering angle cos "& is connected to the invariant t: 

s\ 2 



cos? 



where 



/3/ = \ 1 



4 mj. 



(34) 



(35) 



(36) 
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2.5 The 4-leg CC processes fif[ (W) ff 

In version 1 . 00 we have implemented particular 2/ 2/ CC processes, having in mind their application 
to Drell-Yan type CC processes at hadron colliders as well as for 3-particle top decays. 



u + — > / + Ul 



I e 



Xw{s) 



4s 



7m (1 + 75) «) 7m (1 + 75) ^ll{s, t, u) 



+7m (1 + 75) «) (1 + 75) {-iD^)J^^4s, t, u) 
There are only two non-zero HAs for the case when only one mass (lepton) is not neglected: 



(37) 



W+ — + 
^+-++ 

W+++± 
here and below 



p+ 



= -e^ sini?diP+Xw(s) 



nil 



m. 



J^LL{s,t,u) +1 -] J^LD{s,t,u) 



7i++-± = H+-±- = 'H-++± = H-+-± = H — +± = H ± = , 



P+ = P+{s,mi,0) = ^f7^ 



mr 



4 s 



7m (1 + 75) ® 7m (1 + 75) ^LLis, t, u) 



+7m (1+75) O (1 - 75) {-iDfi)^Ru{s,t,u) 



(38) 
(39) 



(40) 



H+ — + 

n+ — 



P+ 

{I + cos i}uu) —j= Xw{s)J^LL{s,t,u) 



= -e^ sint?„^P+Xiy(s) 



mi 



m. 



J^LL{s,t,u) +1 ^ J^RD{s,t,u) 



H++-± = H+-+± = H-++± = 7i-+-± = H — +± = H ± = . (41) 

Note our angle convention: 'd^i and 'duv are chosen to be the angles between particle momenta in 
the initial and final states in the cms reference frame, and here 

1 



^"^^'^ 2sl s-Ml+iM^r^- 



t{p2) ^ b{pi) + 1+{P4) + l^lips) 

For this case there are four different structures and scalar form factors if the mass of the b quark is 
not neglected. The CA reads: 



7M(l + 75)0 7M(l+75)^i.i.(s,i) 

+7m (1 - 75) <^ 7a. (1 + 75) ^Hi(s, t) 

+ (1 + 75) ® 7m (1 + 75) HD^)T^^{s, t) 

+ (1 - 75) O 7m (1 + 75) {-iDu)J'nn{s, t) 



(43) 
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Here and 4-momcntum conservation read: 



D^, = {pl+P2)^l, P2=Pl+P3+P4: 



while the invariants are 
and 

The four non-zero HAs are 
'^++-+ 



-(P3+P4)^ t = -{pi+P4f, 



dw{s) 



1 



1 



2s2 S-MI + iM^T„ ■ 



(44) 
(45) 
(46) 



H + 

H+ — + 
H-+-+ 



+^ dw{s) s.m.'di^ [P^mt-b - P mt+b) J^ll - {P^mt-b + P mt+b) J^rl 
-P+P- [ (P+ + P-) + (P+ - P-) } , 

dw{s) sini?(| [P+mt-b + P~mt+b) Tll - (P+m*-!, - P~mt+b) J^rl 
+P+P- [ (P+ - P-) + (p+ + P-) ^^^] } , 

+^e^d^{s) (l-cost90^{ (P+-P-):^^^- (P+ + P-)^K.}' 
-^e2d^(s)(l + cost?i)Vi{ (P++P-) J-^^- (P+-P-)^«^}. (47) 



Here 



(mt ± mb)2 - s . 



(48) 



Finally, -ffi is the angle between leptonic 4-monientuni in R-frame {pz + P4 = 0) and the z-axis is 
chosen along pi momentum in the rest frame of decaying top. It is related to t invariant by 



2 1 
t = mt + - 



rrif — ml — s — ^ X{s,m'^,ml) cos i?; 



where X{x,y,z) is the ordinary kinematical function 

X{x, y, z) = + y"^ + — 2xy — 2xz — 2yz . 



(49) 



(50) 



t^b + l- +9i 

This case is similar, although not identical to the previous one. For exact expressions see relevant 
module in the SANC tree: 

EW — > Processes — > 4 legs — > Charged Current — > t— > b 1 nu (HA). 
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2.6 Bhabha scattering 

The CA for Bhabha scattering can be derived from Eqs. (29-30) as foUows (if the electron mass is ne- 
glected): 

^Bhabha = A(S) + A^{s) - [A^{t) + A:,{t)] 



+i e 



S IH. IH' ^ 

.2 Xz{s) 



X I (4^^) '7^ (1 + 75) ® 7m (1 + 75) ^^LLis, t, u) + SJPji, ® 7m (1 + 75) ^ql{s, t, u) 
+iPSe"ff, (1 + 75) O 7m ^lq{s, t, u) + (5g7^ ® 7^ J^qq{s, t, u) 



X {(4'^) 7m (1 + 75) ® 7m (1 + 75) s, u) + Sji^hi. ® 7m (1 + 75) ^q4*, s, u) 

+/f )(5e7^ (1 + 75) <^ 7m ^^Q{t^ s, u) + 6lj^ (g) 7^ J^QQit, s,u)\. (51) 



It is described by the electromagnetic running coupling constant and four FFs with exchanged arguments 
s and t. 

There are six non-zero HAs: 



n = -2e^j{j'^^{t,s,u)-x4t)Se[^!Q{t,s,u)-S,J^^^{t,s,u)]}, 

n+-+- = -e^ (1 - cos6») ^^J^^Q{s,t,u) - Xz{s)S^ [j^^Q{s,t,u) - 5^J^^Q{s,t,u) 
n-+-+ = -e^ (1 - cos^) ^^J^^Q{s,t,u) - Xz{s)5^ [j^^^{s,t,u) - 6^J^^Q{s,t,u) 

= -e^ (I + C0S6I) |j^2Q(s,t,u) 

+Xzis) {KAs.t,u) - 5^ [:Fl^{s,t,u)+:F'^^{s,t,u) - 5^T^^^{s,t,u)\) 
+ \ [•^2q(^' + Xz{t) s, u) - 6^ [^f„(t, s, u) + T^,{t, s, u) - (5e^QQ(i, s, «)] ) ] }, 

H-++- = -e^ {l+cos9){j^^^{s,t,u)+Slxz{s)J'^Q{s,t,u) 

+ J [-^q(*. s, u) + 6lxz{t)J'^^{t, s, «)] }, 
H+++- = H++-± = H+-++ =H+ =0, 

n-+++ = H-+—=H—+± = H—+ = Q, (52) 

however, since for Bhabha scattering = Tqj^-, the number of independent HAs is actually reduced to 
four as expected. 
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2.7 ffbb processes 



In SANG V. 1 .00 we have implemented three classes of ffbb processes: ffjj 0, ffZj and 
f fH^ — > 0. Due to space shortage, we limit ourselves in this paper to the process ffHj 0. Moreover, 
for ffbb processes the variety of cross channels is more rich than for the case of ffff 0. For 
example, for process f fHj ^ it is worth considering at least three channels: 

• annihilation, // H'j; 

• decay, H ffj; 

• and H production at je colliders, — > He. 

For all the channels we still can write down almost unique CA. Below we give it exactly for the 
annihilation channel, f{p2)f{pi) ^ H{p4)"f{p3), but it might be easily rewritten into any other channel. 
This is not the case, however, for the HAs. The latter must be recomputed for any given channel. 

There are eight structures transverse in photonic 4-momentum, 4 vector and 4 axial ones 



9 ^jp- \v{Pi) « (P2) + m)){p2\ - + m)){p,),] eZ{ps)F,,{s, t) 



V {Pl)i)3luU {P2) eZ{P3)Fv2is, i) 

1 



-v{p{) i 

-v{pi) i 



1 



u{p2)£l{Pi)Fv3{s,t) 

u{p2)eZ{P3)Fvi{s,t) 



+v{pi)j!iu{p2) [{U^+m}){p2). - {T^+m}){pi)^] eZ{p3)Fai{s,t) 



-V {piihlulbU {j>2) sZ{j>i)Fa2{s, t) 



-v{pi) i 
-V (pi) i 



"f5U {p2) eZ{p3)Fa3{s,t) 
l5U{p2)eZ{P3)FaA{s,t) \, 



each multiplied by the corresponding FF: F^i-=-^4 and -Fai-=-a4- In above expressions 

T^ = {P2+Psf, U^ = {P2+PAf. 

The eight different HAs are 



H 



+++ 



sPf {l3fF,i{s,t) - Fai(s,t)) - 13+ (F,2(s,t) - Fa2{s,t)) 



+mf {F^3{s,t) - PfFa3{s,t))+mf {F^4{s,t) + PfFa4{s,t)) 



n. 



= ko 



Spf {PfF^i{s,t) + Fai{s,t)) - (3+ {F^2{s,t) + Fa2{s,t)) 



+mf {Fy3{s,t)+/3fFa3{s,t))+mf (F„4(s,t) - l3fFai{s,t)) 



H 



++- 



= fcn 



sPf {PfF,i{s, t) - Fal{s, t)) - 13^ {F„2{S, t) + Fa2{s, t)) 



+mf {Fy3{s,t) - pfFa3{s,t))+mf {Fy4{s,t) - pfFa4{s,t)) 



(53) 



(54) 
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■++ 



sPf {l3fF,^{s,t)+Fai{s,t)) - 13. {F,2{s,t) - F,2(s,t)) 



+m/ (F,3 (s, t) + l3fFa3{s, t))+mf (F„4(s, t) - PfF^iis, t)) 



4^ {F.2{s,t) - F,2{s,t)) 



-P+ {F^3{s,t) + PfFa3{s,t)) - /?_ {F^4is,t) + PfFa4is,t)) 



n-+- = -k^ 



4^(F„2(.S,t)+Fa2(s,t)) 



4^(F„2(s,t)+-Fa2(s,t)) 



-/?_ (j;3(s,i) +/3/Fa3(s,i)) - /3+ (F^4(s,t) +/3/J^a4(s,t)) 



4^(F,2(s,i)-^^a2(.S,i)) 



-/?_ (F,3(s,f) - /3/F„3(s,i)) - /3+ (F,4(s,t) - /3/F„4(s,t)) 



(55) 



with the coefBcients 



ko = -eg 



k± 



Furthermore, /3/ as in Eq. (36) and 



-eg 



Q^1±C0S^^^_^2)^_ 



4\/2 



/3± = l±/3/. 



(56) 



(57) 



The angle is the cms angle of the produced photon (angle between p2 and ps). 

For the sake of completeness we also present the amplitude in the Born approximation. In terms of 
structures (53) it reads: 



A Born 



-eg 



Qf^f I 

{T^+m})iU^+m}) 



x^;(pi)U(C/'+m^)(p2).-(T2 + m^)(pi).] +^(Q' + M2)^37,U(p2)£:^(P3), (58) 



where 



= (Pl + P2f = -S. 



(59) 
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3 Precomputation 



3.1 Introduction 



Ci^ Root 
t lE. SANC 

'h'S QED 

If 13 Precomputation 



Self 

■3 Photon 
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This section is devoted to a rather detailed description of precomputation in SANC. The concept of pre- 
computation is very important for the SANC project (see Ref. [2]). The basic idea here is to precompute as 
many one-loop diagrams and derived quantities (like renormalization constants, various building blocks 
etc.) as possible since the CPU time needed is in general quite large for the above quantities making it 
impractical to compute them in each SANC run. 
Recall our particle notation conventions: 

— / stands for any fermion (lepton or quark); 

— b stands for neutral bosons A, Z, H; 

— when we need to be more concrete, we use I for leptons instead of /, and precisely A, Z, W, H for 
bosons. 

It is worth emphasizing that at the precomputa- 
tion phase it is not necessary to distinguish the pro- 
cess channel. While computing one-loop diagrams 
all 4-momenta are considered as incoming. In the de- 
rived expressions (say for the scalar form factors) any 
required channel is obtained by means of an appro- 
priate permutation of arguments (say of Mandelstam 
variables s, i, u). 

The Fig. 2 shows the fully open menu for "Pre- 
computation" in the QED branch of SANC. 

It consists of Self (Energies), Vertex and Box 
submenus. Self energies, in turn, are subdivided into 
Photon and Lepton submenus. They are further 
subdivided into calculation of diagrams themselves: 
Photon Self and Lepton Self. Precomputed and 
stored one-loop diagrams are used for calculation of 
corresponding renormalization constants: CalcPhot- 
RenConst and CalcLepRenConsts, which are also 
stored. Altogether they arc used for the calculation 
of rcnormalized propagators: CalcPhotRenProp 
and CalcLepRenProp. The latter is used to cal- 
culate the rcnormalized self energy 4-leg diagram for 
Compton scattering in QED UAA RenSelf. The 
file llbb Bornlikect computes Born-like countert- 
erms (see Section 3.2.3). 

Vertex consists of 3-photon-leg AAA, photon-2 
lepton All and 4-leg vertices for NC UAA Compton- 
like QED processes (any channel, see Section 3.3.3). 

Box is represented by the NC 2-photon exchange 
4-lepton-leg box (direct and crossed) 1111 Box (see 
Section 3.4.1) and by the two topologies (T2 and T4) 
of boxes appearing in Compton-like processes llAA 

Box (see Section 3.4.2). Finally, the file AA FF Box transforms results obtained by 1111 Box into the 
scalar form factors of a 41 process. 



9 Q 



<3 Ren Blocks 

fh^ CalcPhotRenConst 
^ CalcPhotRenProp 
¥h^ Photon Self 
Lepton 

Ren Blocks 

- ^1 CalcLepRenConst 

- CalcLepRenProp 

- ¥m IIAA RenSelf 
^ IIAA Bornlikect 

Fct Lepton Self 
Vertex 

fhn AAA Vertex 

All Vertex 
^ IIAA Vertex 
Box 

Mil Box 
' fh^ IIAA Box 
^ AA FF Box 
Processes 



□ 

□ B/V 

□ QCD 



Figure 2: Precomputation in QED part. 
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The tree of "Precomputation" in the EW part, 

® Root 
9 lE; SANC 

ts- 113 QED 
C3 BA/ 
7 ^3 Precomputation 
? <3 Self 

^ <3 Boson 

t G Ren Blocks 

Fct CalcBosRenConst 
Ryi CalcBosSEfunctlons 
Ryi CalcBosct 
^ ffHZ RenSelf 
Boson Self 
R?^ Tadpole 
Fernnion 
G Ren Blocks 

Frs^ CalcFerRenConst 

CalcFerRenProp 
Fi^ ffbb RenSelf 
F tm ffbb Bornllkect 
Fri^ Fermion Self 
I f- Vertex 
■a bbb 
9 G Boson 

F ts^ bbb Vertex 
L ^ bWW Vertex 
^3 Fermion 

R?i bbb Vertex 
□ ffbb bbb 
If- bff 

Fct bff Vertex 
Rm ffbb Vertex 
1^ ffHZ bff Vertex 
If- ^3 Box 
t 3 ffff 

^3 Neutral Current 
NC ZA Box 
f ^ NC ZZ Box 
' fh^ NC WW Box 
NC FF Box 
^3 Charged Current 
CC INI Box 
Fi^ CC FIN Box 
if-'ia ffbb 

9 Q Neutral Current 
^ NC T1 T3 Box 
NC T2T4 Box 
NC T5T6 Box 
F^ NC T7 Box 
Fh^ NC Cont Box 
S Charge Current 
o- Gl bbbb 
113 Processes 
QCD 



which is shown in Fig. 3, has many more branches. 

It also consists of Self, Vertex and Box 
sub-items. 

Self energies are subdivided into Boson and 
Fermion submenus. They are further subdi- 
vided into calculation of diagrams themselves: 
Boson Self and Tadpole (section 3.2.1) and 
Fermion Self (section 3.2.2). Precomputed 
and stored one-loop diagrams are used to cal- 
culate the corresponding renormalization con- 
stants: CalcBosRenConst and CalcFerRen- 
Consts, which are also stored. They are all 
used to calculate the ingredients of renormal- 
ized bosonic propagators: CalcBosSEfunc- 
tions and CalcBosct and renormalized fermi- 
onic propagators CalcFerRenProp. The lat- 
ter is used to calculate the renormalized self en- 
ergy 4-leg diagram for ffbb processes by ffbb 
RenSelf; ffbb Bornlikect computes Born- 
like counterterms (section 3.2.3). 

Vertex consists of 3-boson-leg bbb and 
boson-2-fermion bff vertices. Vertices bbb con- 
tain Bosonic and Fermionic components and 
bbb vertices for ffbb processes (closed here). 
The former is further subdivided into any neu- 
tral leg bbb Vertex and in particular any 
neutral boson to W^W~ bWW Vertex (sec- 
tion 3.3.2). Vertices bff are subdivided into 
any bff 3-leg vertices and 4-leg vertices for 
NC ffbb processes (section 3.3.1). 

Box is subdivided into three large classes: 
ffff, ffbb and bbbb each of them is subdi- 
vided further into NC and CC boxes. The ffff 
class contains a rich collection of 4-fermion-leg 
NC and CC boxes, direct and crossed, (sec- 
tion 3.4.1). The ffbb family is presented by 
now by seven topologies of boxes (NC Tl— 
T7 Box) appearing in NC ffbb processes (sec- 
tion 3.4.2). 

The file NC FF Box transforms results 
obtained by ffff into the scalar form factors 
of a 4/ process, while NC Cont Box reali- 
zes some further manipulations with NC f fbb 
boxes (section 3.4.2). 

The files intended for Charged Current 
boxes for ffbb processes and for bbbb boxes 
exist but are not added for the time being. 



Figure 3: Precomputation in EW part. 
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3.2 Self energies 
3.2.1 Bosonic self energy 

Self energies are the simplest one-loop diagrams. There are three topologies of bosonic self-energy: 

Q q 





-q-pi -q-pi 

i) Bosonic component ii) Fermionic component 

Figure 4: Two point diagrams. 

q 




Pi r- - -p^ 

iii) Bubbles 
Figure 5: One point bosonic diagrams. 

q q 




i wwwxaAaaaaaaA/ i 




Pi 



-Pl 



i ^AAAAAAAAAAAAAAAy J 



Pl 



V 



-Pl 



iv) Bosonic component v) Fermionic component 

Figure 6: Tadpoles. 

For the first and the third topologies it is useful to distinguish bosonic and fermionic components. 

All these diagrams are precomputed and stored in the BS.sav file. At user request they may be 
recomputed by a FORM program accessible via the menu sequence EW — !■ Precomputation Self 
— > Boson Boson Self. The five types of self energies are defined in the specific procedure Diagb(i,j) 
(see the definitions of procedure types in the Introduction). 

To give an example of what the precomputation really docs, we present a simplified version of a very simi- 
lar module accessible via the menu sequence QED Precomputation Self Photon Photon 
Self. It starts with the defining expression vacpol for the photonic vacuum polarization, Fig. 4 ii) and con- 
tinues in five steps of calls to the intrinsic procedures (i) FeynmanRules, (ii) Diracizing, (iii) GammaTrace, 
(iv) Reduction and (v) Scalarizing described in Section 4.2. After steps 1-4) we show the intermediate 
results and after step 5) — the final result. 
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#include Declar.h 
#call GlobalsO 

g vacpol = -Tr*vert (+1 , 12, -12,inu, ii) *pr (12 ,q, ii) 

*vert(-l,12,-12,nu,ii)*pr(12,q+pl,ii)*int; 

(i) 

#call FeynmainRules (0) 
vacpol = 

+ den(l,niel,q)*deii(l,niel,q+pl)*e~2*qel~2*int*Tr* 

(-gd(ii ,mu)*gd(ii ,q)*gd(ii ,nu) *i_*mel-gd(ii ,mu) *gd(ii ,q) *gd(ii ,nu) *gd(ii ,q) 
-gd(ii ,mu) *gd(ii ,q) *gd(ii ,nu) *gd(ii,pl)+gd(ii ,mu) *gd(ii ,nu) *mel~2 
-gd(ii ,mu) *gd(ii ,nu) *gd(ii ,q) *i_*iiiel-gd(ii ,inu) *gd(ii ,nu) *gd(ii ,pl) *i_*mel ) ; 

(ii) 

#call Diracizing(O) 
vacpol = 

+ den(l,mel,q)*den(l,mel,q+pl)*e~2*qel~2*int*Tr* 

(-2*gd(ii ,mu) *gd(ii ,all) *qf (all) *qf (nu)-gd(ii ,mu) *gd(ii , all) *gd(ii ,nu) *qf (all) *i_*mel 
-gd(ii ,inu) *gd(ii , all) *gd(ii ,nu) *gd(ii ,pl) *qf (all)+gd(ii ,mu) *gd(ii ,nu)*inel~2 
+gd(ii ,mu) *gd(ii ,nu) *q.q-gd(ii ,mu) *gd(ii ,nu) *gd(ii ,pl)*i_*mel 
-gd(ii ,mu) *gd(ii ,nu) *gd(ii ,all)*qf (all) *i_*inel ) ; 

(iii) 

#call GcunmaTrace (1) 

vacpol = 

+ den(l ,inel,q)*den(l,mel,q+pl)*e~2*qel~2*int* 
( 4*d_ (mu,nu)*mel"2+4*d_(mu,nu) *q.q+4*d_ (mu.nu) *q.pl 
-4*qf (mu)*pl(nu)-8*qf (inu)*qf (nu)-4*qf (nu)*pl(inu) ) ; 
(iv) 

#call Reduction (0) 
vacpol = 

+ pl(mu)*pl(nu)*e"2*qel"2* 

(-8*b21(pls,l,mel,l,mel)-8*bl(pls,l,inel,l,inel)) 
+d_ (mu,nu)*e~2*qel~2* 

(4*a0 (1 ,inel) -8*b22 (pis , 1 ,inel , 1 ,mel)+4*bl (pis , 1 ,mel , 1 ,mel) *pls) ; 
(v) 

#call Scalarizing(O) 
vacpol = 

+ pl(mu)*pl(nu)*e~2*qel~2* 

(-4/9-8/3*pls"-l*mel~2-8/3*a0(l,inel)*pls~-l+4/3*b0 (pis , 1 ,inel, 1 ,mel) 
-8/3*b0(pls, l,mel, l,mel)*pls"-l*inel"2 ) 
+ d_(mu,nu)*e~2*qel~2 * 

( 4/9*pls + 8/3*inel"2 + 8/3*a0(l ,mel) - 4/3*b0(pls , 1 ,mel, 1 ,mel)*pls 
+8/3*b0(pls,l,mel,l,mel)*mel~2 ) ; 
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We continue describing the structure of the basic program Boson Self: 

The calculation of all bosonic self-energies (diagram-by-diagram) is done by two calls to specific procedure 
CalcBos(i,j), with {i = 3, 6, j = i}, corresponding to the charged external bosons and four calls 
to CalcBos(i,j) with {i = 1,2,4,5, j = i}, corresponding to 7, Z, H and and one call to CalcBos(i,j) 
with {i = 2, j — 1}, corresponding to external bosons Z and 7, respectively. 

Procedure CalcBos(i,j) calls specific procedure Diagb(i,j) and seven intrinsic procedures: 
(i) FeynmanRules, (ii) GammaRight, (iii) Diracizing, (iv) GammaTrace, (v) Reduction, (vi) Sing and (vii) 
Scalarizing. The intrinsic procedures are described in Section 4.2. 

At each call to procedure CalcBos(i,j), the diagrams Figs. 4, 5, 6 are calculated for twelve virtual 
fermions field indices k, I running over 11, 12 ... 22 (corresponding to Ve, e, u, d, ly^, /x, c, s, Vr, t, t and 
h, respectively) and for ten virtual bosons k.l = 1,2, . . . , 10, corresponding to 7, Z, W^, H, (fp, (f)"^ and 
four Faddeev-Popov ghosts X~^,X~, Yz, Ya, respectively. 

The topology of the self-energy diagrams is specified in procedure Diagb(i,j) in terms of the vertices 
and propagators of the diagrams. In the class of diagrams available at present, vertices are of two kinds: 
boson-fermion-fermion (Bff) and three-boson (EBB) vertices. The diagrams are computed in nested loops 
over all allowed field indices of the virtual particles. 

Two more FORM programs are accessible via menu sequences 
QED — ^ Precomputation Self Photon Boson Self — > Photon Self and 
EW Precomputation Self Boson Tadpoles. 

They have very similar structures and compute respectively photonic self energy (vacuum polarization) 
in the QED tree of SANG and tadpole diagrams separately from self energies. The results are stored in 
PSqed . sav and TP . sav files, respectively. 

Precomputed bosonic self-energies are used by the other FORM programs which calculate bosonic 

counterterms 

EW Precomputation — ^ Self Boson — > Ren Blocks CalcBosRenConst 

and the photonic counterterm 

QED — > Precomputation — > Self — > Photon — > CalcPhotRenConst 

in the QED tree; bosonic self energy functions 

EW Precomputation Self Boson Ren Blocks CalcBosSEFunctions, 

and the renormalized photonic propagator in the QED branch 

QED Precomputation — > Self — > Photon CalcPhotRenProp. 

Their results, in turn, are used by a FORM program which computes counter term blocks (crosses), see 
Fig. 15 of Ref. [9]. The latter is accessible via the chain 

EW — > Precomputation — > Self — > Boson — > Ren Blocks — > CalcBosct. 
3.2.2 Fermionic self energy 

There are two topologies of fermion self energy diagrams (two point diagrams Fig. 7 and tadpoles Fig. 8). 



q 




Pi 



~q-pi 



-Pi 



Figure 7: Two point diagrams. 
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ii) Bosonic component iii) Fermionic component 



Figure 8: Tadpoles. 

These diagrams are precomputed and stored in the FS . sav file. They may be recomputed by a FORM 

program accessible via the menu sequence 

EW Precomputation Self Fermion Fermion Self. 

The three types of diagrams are defined in the specific procedure Diagf(i,j). 

Structure of the program Fermion Self: 

The calculation is done by 12 calls to specific procedure CalcFer(i,j), with i = j = 11, 12 . . . , 22, corre- 
sponding to all fermions of four generations, respectively. 

Procedure CalcFer(i,i) calls specific procedure Diagf(i,j) and seven intrinsic procedures (i) Feynman- 
Rules, (ii) GammaRight, (iii) Diracizing, (iv) GammaTrace, (v) Reduction, (vi) Sing and (vn) Scalarizing, 
see Section 4.2. 

At each call to procedure CalcFer(i,i), the diagrams Figs. 7 and 8 are calculated for twelve virtual 
fermions i = 11, 12 . . . , 22 and only six virtual bosons ? = 1, 2, . . . , 6, corresponding to 7, Z, W^, h, (jP 
and (j)"^ , respectively, since Faddeev-Popov ghosts do not contribute here. 

The topology of the self-energy diagrams is specified in procedure Diagf(i,j) in terms of the vertices 
and propagators of the diagrams as well as sign and combinatorial factors (sign('r) and cft(T)). 

There is one more FORM program, accessible via the menu sequence 
QED — > Precomputation — > Self — > Lepton — > Lepton Self 

whose structure is very similar to program Fermion Self. It computes the leptonic self energy in the QED 
tree of SANC. The results are stored in LSqed.sav file. Precomputed fermionic self-energies are used by 
FORM programs which calculate fermionic counterterms 

EW — > Precomputation Self Fermion Ren Blocks CalcFerRenConst 

and leptonic counterterms 

QED Precomputation Self — > Lepton — ^ CalcLepRenConst. 

Both fermionic self-energy diagrams and fermionic counterterms are used by FORM programs which 

compute the renormalized fermionic (leptonic in the QED branch) propagators 

EW Precomputation Self Fermion Ren Blocks CalcFerRenProp and 

QED Precomputation — > Self Lepton CalcLepRenProp, respectively. 

In the FORM code CalcLepRenProp we show the main steps of the calculations. The two latter 
codes end up demonstrating the vanishing of the renormalized propagator on the corresponding fermionic 
mass shell, as is required by the on-mass-shell (OMS) renormalization scheme, see Section 2.2 of [9]. 
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3.2.3 Fermionic self energy for ffbb processes 



Moving to precomputation of the building blocks for f fbb processes, we change our conventions. Now any 
object: self energy, vertex and box are considered to be 4-legs, rather than 2-, 3- and 4-legs respectively, 
as we did before. The main motivation for this change is our observation that vertices with off-shell 
fermions are inconvenient to treat and the resulting expressions are more compact if we consider 4-legs 
on-mass-shell building blocks instead of 3-legs off-shell. Accepting this convention for vertices, it is natural 
to treat self-energies also like 4-leg objects shown in Fig. 9: 

fu=typcIU fu=typeIU 



Here the blob at the fermion propagator denotes the sum of all self-energy diagrams described in 
Section 3.2.2. These self energies are precomputed by a FORM program accessible via the menu sequence 
EW — > Precomputation — > Self — > Fermion — > Ren Blocks ffbb RenSelf. 

The required CPU time is still very short, and at user request they may be re-computed. The two types 
of self energies are defined in the specific procedures CalcFerSEt(fu,fd,vd,vu) and CalcFerSEu(fu,fd,vu,vd), 
which use renormalized fermionic propagators precomputed by CalcFerRenProp. Note that the labels 
't' and 'u' are associated with the Mandclstam variables t,u, see Section 2.4. 

Both t and u channel procedures call intrinsic procedures FeynmanRules, GammaRight and p2l where 
the latter, together with several id's in between, performs obvious identities and change of variables. 

Computed diagrams FSEt'fu"fd"vd"vu"jl"kl' and FSEu'fu"fd"vu"vd' jl"kr are stored in the file 
ffbbSelfxi'xi"fu"fd"vu"vd'.sav, where the predefined parameter 'xi' has the following meaning: 
#define xi "0" / * .eq.O to test gauge invariaince / * .eq.l to work in xi=l gauge 
This option is introduced mostly to save CPU time since calculation in the xi=l gauge are much faster 
than in the gauge. However, it is always appealing to see the explicit cancellation of gauge parameters. 
That is why we try to maintain option xi " 0" as long as possible even though in some cases it is extremely 
time consuming. 

There is a similar FORM program in the QED part, accessible via the menu sequence 
QED Precomputation Self Lepton Ren Blocks llAA RenSelf. 

Furthermore, the code, accessible via menu sequence: 
EW — !■ Precomputation Self — ^ Fermion Ren Blocks ffbb Bornlikect 
computes contributions from four diagrams similar to Fig. 9, in which the self energy blob at the fermion 
propagator is replaced by counterterm 'crosses' (one for each of the four diagrams at each vertex with 
indices fi and v, similar to Fig. 12). The result is stored in the file fFbbBornlikectxi'xi"fu"fd"vu"vd'.sav. 

The file from QED Precomputation Self Lepton Ren Blocks ffAA Bornlikect 
does the same job in the QED part. 
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Figure 9: Self energy ffbb diagrams. 
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3.3 One-loop vertices 



The current SANC version has all bff and bbb 3-leg SM vertices. 



3.3.1 bff vertices 

There are two topologies of bff vertices: FBF and BFB: 
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Figure 10: bf f vertices. 



These vertex diagrams are precomputed by a FORM program accessible via menu sequence 
EW Precomputation Vertex bff — * bff Vertex. 

At user request they may be recomputed. The two types of vertices are defined in the specific 
procedures FBF(iJ,l) and BFB(i,j,l), where the topologies of all vertices of the type of Fig. 10 are specified 
in nested loops over all allowed field indices of the virtual particles. 
Structure of program bff Vertex 

The calculation is done by a single call to specific procedure CalcVertex('typeB','typeU','typeD'), followed 
by a call to intrinsic procedure p2D. 

The arguments typeB, typeU, and typeD are the field indices of the external particles, set in the 
command that starts this program: typeB is the incoming boson, typeU is the antifermion with incoming 
4-momentum p2, and typeD is the fermion with incoming 4-momentum p\. For instance, to calculate the 
vertex for Z ^ uu decay, typeB, typeU, and typeD are set equal to 2, 13, 13, respectively. The results are 
stored in files V typeB' 'typeU' 'typeD' . sav. 

Procedure CalcVertex(i,j,l) calls specific procedures FBF(i,j,l) and BFB(i,j,l) and the intrinsic procedures 
in the following sequence: FeynmanRules, GammaRight, Diracizing^, Reduction, Diracizing, Pulling, Diraceq, 
Sing, Scalarizing, 2Qs, Masshell and Symmetrize. The intrinsic procedures are described in Section 4.2. 

In the QED tree there is a similar entry 

QED Precomputation Vertex Aff Aff Vertex, 

which computes 7/! vertices for a lepton of a kind j=typeU=12, 16,20. These vertices are defined in the 
specific procedure DiagVert(j). The results axe stored in files Vqed'typeU' .sav. 

^For W ff vertices, this procedure is called exceptionally in the so-called special use mode, see Section 4.2. This is due 
to a delicate treatment of the auxiliary Passarino-Veltman functions originating from a double pole l/(p^)^ in the photonic 
propagator in the gauge. 
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3.3.2 bbb vertices 



The bosonic component of three-boson vertices has four topologies shown in the following diagrams: 
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Figure 11: Four topologies for three-boson vertices. 
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These vertex diagrams are precomputed by a FORM program accessible through the chain of chcks 
EW — > Precomputation — > Vertex — > bbb — > Boson bvv Vertex. 

They also may be recomputed at user request. The four topologies are defined in the specific procedures 

TribosVertTl(i,l,j). 

Structure of program bvv Vertex 

The calculation is done by a single call to specific procedure CalcTribosVertT14('typeB','typeD','typeU') 
which calls first of all four topologies in Figs. 11 defined in specific procedures TribosVertTk(i,l,j), with 
k = 1,2,3,4. 

Just after that four specific procedures CITribosVertTk(i,IJ) perform dusterizing of the computed 
diagrams. One should note that clusters have different meanings in SANC. Here dusterizing means nothing 
but summation over all virtual field indices, resulting in the dependence of cluster names only on external 
field indices. 

After dusterizing, the usual calls of the intrinsic procedures follow: FeynmanRules, 2Qs, Reduction, 
2Qs, Sing and Scalarizing. 

The fermionic component of three-boson vertices has only one topology Tl. The diagrams are defined 
in the specific procedure TribosVertf(i,l,j). 
Structure of program Fermion — > bbb Vertex 

The calculation is done by a single call to specific procedure CalcTribosVert('typeB','typeD','typeU') which 
at first calls the topology defined in the specific procedure TribosVertf(i,l,j). Next, specific procedure 
CITlfer(i,l,j) performs dusterizing. 

The calculation of the cluster is followed by calls of the intrinsic procedures: FeynmanRules, Gannma- 
Left, Diracizing, GammaTrace, Reduction, Scalarizing, Sing, Scalprod and DivisionGramDet. 

One may access a similar FORM code to precompute the HWW, ZWW and AWW vertices via 
sequence EW — > Precomputation Vertex — > bbb — > Boson bWW Vertex. 

3.3.3 Vertices for ffbb processes 

There are four blocks of vertices met in ffbb processes where only diagrams with fermion exchange 
contribute at the Born level (we have also ffZH process where there is the Born diagram with boson 
exchange, so-called Higgsstrahlung, but this process is not added to SANC v. 1.00): 

The four building blocks are precomputed by a FORM program accessible via menu sequence: 

EW Precomputation Vertex — > bff ffbb Vertex. 

Structure of program fFbb Vertex 

The calculation is done by 2 (g) 3 calls to the following specific procedures: 
CalcVert('typeFU','typelD','typelU',mu) — computes vertex blob with index 

CalcVertmut('typelU','typelD','typeFD','typeFU',mu) computes diagram b); 
CalcVertmuu('typelU','typelD','typeFU','typeFD',mu) — computes diagram c); 
CalcVert('typeFD','typelD','typelU',nu) — computes vertex blob with index v; 

CalcVertnut('typelU','typelD','typeFD','typeFU',nu) — computes diagram d); 
CalcVertnuu('typelU','typelD','typeFU','typeFD',nu) — computes diagram a). 

The procedure CalcVert recomputes bff vertices for given set 'typeFU'.'typelD'.'typelU' and creates 
expressions vfbf'vu"fd"fu"mu"kl"k2"k3' and vbfb'vu"fd"fu"mu"kl"k2"k3' similar to those described in 
Section 3.3.1, but not applying intrinsic procedure Masshell since in the case under consideration one of 
fermions is off-shell. The four procedures CalcVertmut, CalcVertmuu, CalcVertnut and CalcVertnuu have 
very similar structures. 

In each of four procedures the intrinsic procedures FeynmanRules, GammaRight, bpldentities, p2m, 
p2l and p2p are called. Next follows an important shift of the virtual fermion 4-momenta to 'real' four 
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Figure 12: Four f fhh vertex diagrams. 



momenta pi,p2,P) <3) and finally ExtMomentumWI is applied. 

All the building blocks are clusterized by summation over virtual momenta in vertex loops (see dis- 
cussion of clusterization in Section 3.4.2.). 

Besides xi there are three more internally defined options: on (see Section 3.4.2 for more discussion 
about option on) and mf, mp, of which the latter two allow to neglect fermion masses. 

#define on "1" 

* .eq.O photons are off mass-shell; .eq.l photons are on mass-shell 
#define mf "0" 

* .eq.O fermion mass, pm( 'typelD' )=0; .eq.l it is not neglected 
#define mp "0" 

* .eq.O partner mass, pm('typeIDp')=0; .eq.l it is not neglected 

The result of this calculation is stored in the file 
ffbbVertxi'xi'on'on'mf mf mp'mp' 'typelU"typelD"typeFU"typeFD'.sav. 



27 



3.4 One-loop boxes 

Approaching the description of precomputation of boxes, one should note that the world of boxes is much 
more rich and complex compared to self-energies and vertices. If for the latter case we still could profess 
the idea of allowing recomputation of all one-loop vertices needed for the process under consideration, 
then in the case of boxes we must change our strategy if the user wants to carry out some recomputation. 

For boxes, the idea of precomputation becomes vitally important for realization of SANC project. As 
will be explained below, calculation of some boxes for some particular processes takes so much time that an 
external user should refrain from repeating precomputation. Furthermore, the richness of boxes requires 
more classification. Depending on the type of external lines (fermion or boson), we will distinguish three 
large classes of boxes: ffff, ffbb and bbbb. 



3.4.1 Boxes for //// processes 

The ffff boxes, which are met in the description of 4f processes, are of two topologies, direct and crossed, 
and arc characterized by two fcrmionic currents ii and jj, coupled by two bosons with field indices kl and 
k3; see Fig. 13 showing all the field indices (for external and virtual fields), types of external fermions 
typeXX, Lorenz indices and momentum flows for the direct and crossed topologies. 
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Figure 13: Direct and crossed boxes for //// processes. 

These diagrams describe both neutral current (NC) and charged current (CC) boxes both for 1^ 3 
decays and 2— > 2 reactions in any channel {s,t,u). In the general case, virtual boson field indices run 
from 1 to 6, and virtual fermion indices run over doublets {typelD,typeIDp} and {typeFU.typeFUp}, where 
typelDp and typeFUp are isospin partners of typelD and typeFU, respectively. 

The calculation of direct and crossed boxes is realized in two specific procedures direct(iu,id,fu,fd) and 
crossed(iu,id,fu,fd) in nested loops over all allowed field indices of the virtual particles. 

Such a realization may apparently take into account all the external fermion masses. However, in 
practical applications we usually treat fermion masses of one of two currents, ii or jj, as being massless, 
which is true for the processes of the kind fi + fi + f + f ^ 0, upon which we concentrate. Here fi 
denotes a massless first generation fermion or any neutrino. 

Various FORM codes computing these boxes are accessible via chains of clicks 



28 



EW ^ Box ^ ffff ^ Neutral Current -> NC ZZ(ZA,WW) Box and 
EW ^ Box ^ ffff ^ Charged Current ^ CC INI(FIN) Box. 

Here INI(FIN) means which of the masses, initial fermions ii or final fermions jj are kept non-zero. All 
the codes have very similar structures with minor modifications. The calculation is done by a single 
call to specific procedure CalcBoxNC('typelU','typelD','typeFU','typeFD') which calls first the two specific 
procedures — direct(iu,id,fu,fd) and crossed(iu,id,fu,fd). In those procedures all the box diagrams are 
defined in terms of vertices, propagators and external spinors tlo(pl), tro(p2), tle(p3) and tre(p4). 
The subsequent calculations are very similar for all 4f-boxes, both NC and CC. 

In the beginning of each NC(CC) XX Box program, described in this Section, there is an internal 
definition #define xi "0" which chooses among two internal options with an obvious meaning described 
in Section 3.2.3. 

Note, that for 4f boxes we gain only little CPU time choosing xi "1" option, since action of BoxPre- 
reductionNC(CC) results in nearly complete cancellation of ^ already before Scalarizing. For example, the 
program NC ZA Box needs about 3 minutes CPU time at a 1.6 GHz computer for both options. A similar 
picture is valid also for the program CC FIN Box. 

Actually, 4f-boxes are the last precomputation programs which still may be recomputed at user 
request. However, we do not recommended to recompute them, because 3 min is already a noticeable 
time. 

The results of calculations of NC boxes are stored in the files ffffZZxi'xi"typelU"typelD"typeFU"typeFD' 
.sav and ffffZAxi'xi"typelU"typelD"typeFU"typeFD'.sav and are loaded by a FORM program 

EW Box ^ ffff ^ Neutral Current NC FF Box 

which constructs box form factors stored in the files ffffNCFF'typelU"typelD"typeFU"typeFD'.sav. This 
special FORM program will be described elsewhere. 

3.4.2 Boxes for ffbb processes 

There are seven topologies of boxes which are met in the description of 2f2b processes. Their enume- 
ration is borrowed from Ref. [1] . All these boxes have apparently only one fermionic current conventionally 
marked by the current index ii but different numbers of internal bosonic lines. 

Topologies T2, T4 

Wc begin with the simplest case of topologies having only one virtual bosonic line, sec Fig. 14 showing all 
the field, type and Lorentz indices and momentum flows. These two topologies are actually of the type 
of direct and crossed ones considered in the previous section. However, their structure is quite different. 

These diagrams also describe both NC and CC boxes in any channel (s,t,u), both decays and reac- 
tions. The virtual boson fleld indices run from 1 to 6, and virtual fermion indices run over doublets 
{typelD.typelDp}, where typelDp is the isospin partner of typelD. It is foreseen to cover CC processes in 
the future; currently wc have only NC processes. 

The calculation of T2 and T4 topologies is realized in two specific procedures boxT2(fu,fd,vd,vu) and 
boxT4(fu,fd,vu,vd) in nested loops over all allowed field indices of the virtual particles. Contrary to the 
4f-casc, we do take into account the external fermion masses. 

The calculation starts by two calls to specific procedures CalcBoxT2('typelU','typelD','typeFD','typeFU') 
and CalcBoxT4('typelU','typelD','typeFU','typeFD') which call two specific procedures shown above. 

Then, for each topology, the calculation continues by clustering boxes calling the procedure Clus- 
terboxT2(4)('fu','fd','vd','vu') which creates four clusters, CI'k'T2(4)'fu"fd"vd"vu' for k=l,4, depending 
on ^z, and independent of any ^, respectively. Next follow calls to the intrinsic procedures 
FeynmanRules, GammaRight, Diracizing, Diraceq, Reduction, Pulling, DIraceq, PulllngOrder, Scalprod, Sing, 
ExtlVlomentumWI described in Section 4.2. 
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Figure 14: Boxes for ffbb processes, topologies T2 and T4. 

Then follows a do-loop over k making separate scalarizations of four clusters inside which six more 
intrinsic procedures are called: Scalarizingdp, Scalarizing, DivisionGramDet, bpldentities, p2m, p2p. 

In the beginning of the program there arc four usual definitions: 
#define xi "0" / #define on "0" / #define mf "1" / #define mp "1". 
the meaning of three of which was explained in the previous sections. 

Note that for the 2/2& processes where one (or two) bosons are photons one may, of course, choose 
on "1" which greatly saves CPU time (see tables of CPU time below in this section). Option on "0" is 
basically foreseen for 2^3 processes where these 2/26 boxes with off shell boson(s) are building blocks. 

The results of calculations of these boxes arc stored in the files 
ffbbT2xi'xi'on'on'mf'mf 'mp'mp"fu"fd"vd"vu'.sav and ffbbT4xi'xi'on'on'mf'mf 'mp'mp"fu"fd"vu"vd'.sav. 

The calculation of 2/26 boxes takes a lot of CPU time. The NC T2T4 Box is the fastest, but even 
this program should not be recomputed by users. The other //66 box topologies, except T7, take much 
more CPU time. 

In conclusion of this section we present an instructive example of how much CPU time is needed to 
compute T2+T4 topologies for the dd 77 process at a 3 GHz PC running Linux: 

xi " 0" , on " 0" , mf " 1" , mp " 0" — 5 hours, 

xi " 0" , on " 1" , mf " 1" , mp " 1" — 70 minutes, 

xi " 1" , on " 1" , mf " 1" , mp " 1" — 60 minutes. 

This is why a recomputation of these boxes is not allowed. 

Topologies Tl, T3 

We jump now to the most complex case of topologies having three internal bosonic lines, see Fig. 15 
showing all the field, type and Lorentz indices and momentum flows. These two topologies are also of the 
type of direct and crossed ones, however their structure is much more complex than anything considered 
so far. 

The paragraph after Fig. 14 applies also in this case. 

The calculation of Tl and T3 topologies is realized in two specific procedures boxTl(fu,fd,vd,vu) and 
boxT3(fu,fd,vu,vd) in nested loops over all allowed field indices of the virtual particles. As before, we take 
into account the external fermion masses. 

So far we have considered 2/26 processes for the case when one of the bosons is a photon. In such a 
case the internal bosons must be charged (with field indices ±3 and ±6). 
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Figure 15: Boxes for ffbb processes, topologies Tl and T3. 

The calculation starts by two calls to specific procedures CalcBoxTl('typelU','typelD','typeFD','typeFU') 
and CalcBoxT3('typelU','typelD','typeFU','typeFD') which call the two specific procedures mentioned above. 

Then, for each topology, the calculation continues by calling the procedure ClusterboxTl(3)(fu,fd,vd,vu) 
which creates in this case only two clusters, Crk'Tl(3)'fu"fd"vd"vu' for k—2,3, where k=2 collects all 
neutral virtual bosons and k=3 all charged bosons. Next follow standard calls to the intrinsic procediires 
FeynmanRules, GammaRight, Diracizing, Diraceq, Reduction, Pulling, Diraceq, PullingOrder, Scalprod, Sing, 
ExtMomentumWI described in Section 4.2. 

Then follows a call to the intrinsic procediuc ScalarizingProj which scalarizcs two clusters, k=2,3, 
splitting them first into as many Dirac-Lorentz structures as it has. This is done in order to avoid 
limitations of FORM v3.0 which cannot handle internal files of length greater than 1.6 Gb or so. This 
might be circumvented by switching to FORM v3. 1, however, we did not manage to switch to this version 
so far. 

Inside itself procedure ScalarizingProj creates many intermediate expressions to which intrinsic pro- 
cedures Scalarizing, DivisionGramDet, p2m are applied. At the end, ScalarizingProj collects these pieces 
together again. 

The intrinsic procedure p2p is called at the end. 

The results of calculations of boxes of Tl and T3 topologies are stored in the files 
ffbb'k'Tlxi'xi'on'on'mf mf'mp'mp"fu"fd"vd"vu'.sav and 
ffbb'k'T3xi'xi'on'on'mf mf mp'mp"fu"fd"vu"vd'.sav. 

The label 'k' stands for two clusters: = 2 with all neutral virtual bosons, and = 3 with all charged 
virtual bosons 'kl','k3','k4'. 

In the beginning of the program there are four usual definitions: 

#define xi "0" / #define on "0" / #define mf "1" / #define mp "1" 

having the same meaning as explained in the previous section. 

A table of CPU times for the dd ^ 77 boxes of topologies T1+T3 looks as follows: 
xi " 0" , on " 0" , mf " 1" , mp " 0" — 90 hours, 
xi " 0" , on " 1" , mf " 1" , mp " 1" — 7 hours, 
xi " 1" , on " 1" , mf " 1" , mp " 1" — 14 minutes. 
Of course, a recomputation of these boxes is not allowed. 



31 



Topologies T5, T6 

Now wc switch to the intermediate case of topologies having two internal bosonic lines, see Fig. 16 showing 
all the field, type and Lorentz indices and momentum flows. These two topologies are not the usual pair 
of direct and crossed ones. Note that they are both drawn as two different direct boxes. 
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Figure 16: Boxes for f fbb processes, topologies T5 and T6. 

The paragraph after Fig. 14 is applicable here too. For these topologies, contrary to all previous 
cases, the leptonic current flows through the diagram: from lower-left (p2) to upper-right corner (p4), 
rather then from lower-left (p2) to upper-left corner (pi). This forces us to introduce the notion of in the 
sense of Reduction, i.e. to perform the calculation of a diagram denoting momentum flows as is suitable 
for Reduction and at the end come back to the Real notation for external momenta and Mandelstam 
invariants. This is why the corresponding procedure is called isoR2Real, see Section 4.2. 

The calculation of T5 and T6 topologies is realized in two specific procedures boxT5(vd,fd,vu,fu) and 
boxT6(vu,fd,vd,fu) in nested loops over all allowed field indices of the virtual particles. 

Even though we consider the processes with a photon, for this case internal bosons can be charged or 
neutral if the photon is coupled to the fermion line. The virtual fcrmion index runs over a doublet. 

The calculation starts by several calls to specific procedure CalcBoxT5('typeFD','typelD','typeFU','typelU', 
k3min,k3max,k4min,k4max) and by a single call to CalcBoxT6('typeFU','typelD','typeFD','typelU',k3min,k3max, 
k4min,k4max). Tlicy, in turn, call the two specific procedures shown above. 

Such an asymmetry is due to the fact that we assume 'typeFD' to be a photon, and therefore the 
virtual bosons in diagram T6 must be charged: k3min=3, k3max=6, k4min=3, k4max=6. Note, that 
k3={k3min,k3max} and k4={k4min,k4max}. 

For diagrams with two virtual bosons the clustering is performed in the following way: 

— Cluster 22, k3={2,5}, k4={2,5} 

— Cluster 33, k3={3,6}, k4={3,6} 

— Cluster 42, k3={4,4}, k4={2,5} 

— Cluster 24, k3={2,5}, k4={4,4} 

— Cluster 44, k3={4,4}, k4={4,4} 

Then the calculation continues for each topology by calls to the intrinsic procedures Feyn man Rules, 
GammaRight, Diracizing, Reduction, Pulling, Diraceq, PullingOrder, Scalprod, Sing, ExtMomentumWI and 
then to ScalarizingProj as described in the previous section. 

The intrinsic procedures p2p and isoR2Real are called at the end. 
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The results of ealculations of boxes of T5 and T6 topologies are stored in the files 
ffbb'k3min"k4min'T5xi'xi'on'on'mp'mp"vd"fd"vu"fu'.sav and 
ffbb'k3min"k4min'T6xi'xi'on'on'mp'mp"vu"fd"vd"fu'.sav. 

In the beginning of the program there are usual definitions: 

#define xi "0" / #define on "0" / #define mf "1" /#define mp "1" 

A table of CPU times for the dd 77 boxes of topologies T5+T6 reads: 
xi " 0" , on " 0" , mf " 1" , mp " 0" — 32 hours, 
xi " 0" , on " 1" , mf " 1" , mp " 1" — 4.5 hours, 
xi " 1" , on " 1" , mf " 1" , mp " 1" — 19 minutes. 
A recomputation of these boxes is not allowed either. 

Topology T7 

Topology T7 also has two internal bosonic lines, see Fig. 17, however it is rather a pinch of topologies Tl 
and T4 (bosonic line with field index k4 is pinched out). 

fu=typeIU vu=typeFU 




Figure 17: Boxes for ffbb processes, topology T7. 

The calculation of the T7 topology is realized in the specific procedure boxT7(fu,fd,vd,vu) in nested 
loops over field indices. If there is only one photon in the final state, the virtual bosons must be charged. 
This is why there is only one cluster in this case. 

This is a vertex-like diagram and for this reason it is much simpler then all the others. 

The calculation starts by a single call to specific procedure CalcBoxT7('typelU','typelD','typeFD','typeFU') 
which calls the procedure boxT7 shown above. 

The calculation continues by calls to the intrinsic procedures FeynmanRules, GammaRight, Diracizing, 
Diraceq, Reduction, Pulling, Diraceq, Scalprod, Sing, ExtMomentumWI, Scalarizing, p2m, DivisionGramDet. 

The results are stored in the file ffbbT7xi'xi'on'on'mp'mp"fu"fd"vd"vu'.sav. 

In the beginning of the program there are four usual definitions: 
#define xi "0" / #define on "0" / #define mf "1" / #define mp "1". 

The topology T7 needs only a few seconds of CPU time. 
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4 SANC Procedures 



4.1 Introduction 

At Level 1 the SANC database contains FORM programs and procedures. All FORM programs are acces- 
sible for the user via a sequence of clicks on SANC tree, when we reach a file and open it. The situation 
is very different for the procedures which are of three kinds: specific, special and intrinsic. 

A procedure is called specific if it is used only in one particular program. Normally, specific procedures 
are included in the corresponding FORM program body and are also open for the user. Moreover, they 
are very easy to read and describe. 

Special procedures are usually iiscxl by only a limited number of FORM programs and, similar to 
specific procedures, they do a job which is relevant for these programs only. However, contrary to specific 
procedures we do not open them for the users, since their content is not so transparent to be read and 
described. It is envisaged to upgrade special procedures to the level of intrinsic procedures in future 
versions of SANC. 

Intrinsic procedures are used by many FORM programs and should be easily used in any new program. 
Their functions are totally determined by the list of their arguments which are of two types: genuine 
arguments (FORM variables), denoted below as AVALUE, and options, usually integer numbers, denoted 
as IVALUE, which are used as switches governing calculation fiow inside the procedure. Sometimes an 
IVALUE (see below) stands for a special use, i.e. a special flow inside an intrinsic procedure. 

4.2 Intrinsic procedures 

Below we give a list of intrinsic procedures which are mostly met in the Precomputation trees of 
SANC. In this list, which is not complete by far, the procedures are presented in alphabetic order, and 
the 'treatment' of their options and arguments is explained. 

a2b : replaces symbol "a" to symbol "b". Possible arguments: Sw Cw, — > Mz, 76 — * 75 ) 77 — * 75, 
Oi (Jj, 5i Sj etc. and vice versa. 

AVALUE = (s wjC^y) , for example. 

bpldentities : applies identities to the so-called auxiliary Passarino-Veltman functions (bpO=6o and 

bpl=6i [1]), attempts to exclude bpO (and bpl). 

IVALUE = (I) 

1=0 — to exclude bpO 

1=1 — to exclude bpO and bpl. 

Diraceq : apphes Dirac equations to expressions preliminary simplified with the aid of Pulling. 



AVALUE = (i,j,k,l) 



i for spinor tlo, Dirac equation: tlo^^i 
j for spinor tro, ^2 tro 

k for spinor tre, tle/3 
i for spinor tie, ^4 tre 



—irrii tlo 
irrij tro 
—imk tie 
imi tre 



where (i,j,k,l) are field indices. 
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Diracizing : expressions of the form qVq and ja^ja are simplified; here q = j ■ q, and F is a string of up 
to five 7 matrices. The final step consists of setting q-q = q^ and 7a 7a = n, where n = 4 — e is the 
dimension of momentum space. 

IVALUE = (I) 

1=0 normal use 

1=1 a special use in Wff vertices. 
DirectProdSumm : performs summation in direct products of 7 matrices such as in well known identity 

luloilule O Iull3lnl6 = '^1016 ® 7a 76; 

this procedure knows 76 identities of such a kind. 
AVALUE = (i,j,k,l) , the same arguments as in Diraceq. 

DivisionGramPet : realizes various possibilities to use the algebra of Gram determinants to simplify raw 

expressions. 

IVALUE = (I) 

1=0 division of det3i (active in all subsequent options) 
1=1 division of det4i for Os,Ts topology, for example T1,T2 
1=2 division of det4i for Og, Us topology, for example T3,T4 
1=3 division of det4i for Tg, Ug topology, for example T5 
1=4 division of det4i for Us,Ts topology, for example T6. 

Expansions : expands Bq and Bq functions for small values of some of its arguments. 
AVALUE,IVALUE = (FI) 
FI field index 

ExpansionPhotMassShell : puts an external bosonic momentum in ffbb processes to the corresponding 
mass shell. For example in the process // -yB for or p| its action means 

IVALUE,AVALUE = (I,J,mp,pGs,pBs) 

1=1, J=l pGs = 0, pBs = 
1=1, J=2 pGs = 0, pBs = -mp^ 
1=2, J=l pGs = -mp^, pBs = 
1=2, 3=2 pGs = -mp^, pBs = -mp^. 

ExtMomentumWI : applies Ward identities for external vector boson momenta, i.e. sets {pi)fj, = and 

AVALUE = (I,mu,J,K,nu) 

f2f : realizes the possibility (in particular cases for certain arguments of PV functions) of replacing 
Bo Ao, B^ B^, Bo hi and vice versa, if concrete arguments of the PV functions allow such 
replacements. 

AVALUE = (bO,aO) , for example. 
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Feyn man Rules : applies Feynman rules for propagators and vertices, see Section 4.3. 

IVALUE = (I) 

1=0 for QED part 

1=1 for EW and QCD parts. 

Gam ma Left : all Dirac matrices 75, 76 = 1 + 75 and 77 = 1 — 75 are moved to the left and the expression 
is simplified using identities (75)^ = 1, 767? = 0, etc. 

GammaRight : the same as in GammaLeft, but all matrices 7^ are moved to the right. 

GammaTrace : the traces of products of 7 matrices are evaluated in n dimensional space. 

Globals : performs global declarations by FORM Tables of particle names, particle masses, electric 
charges, ghost charges, mass ratios, coupling constants, weak isospins, gauge parameters, com- 
binatorial factors etc. 

isoR2Real : realizes the ideology of shifting from the level in the sense of reduction pi to the real 4- 
momentapi, {jpi)input {Pi)output and (Invariants)input ^ (Invariants)output, see item Topologies 
T5, T6. 

AVALUE = (plout,p2out,p3out,p4out,Qsout,Tsout,Usout) 
Only output values appear in the argument list; the input string is assumed to be pl,p2,p3,p4,Qs,Ts,Us. 
m2zero : sets a mass mp to in expressions and in the arguments of all functions. 

AVALUE = (mp) 

Masshell : This procedure has four arguments, which must be fermionic field indices; a field, whose index 
is an argument of Masshell, is put on its mass shell. Thus the command #call Masshell ( ' iu ' , , , ) 
sets p('iu')"2 equal to -pin('iu')"2. 

AVALUE = (iu,id,fu,fd) , the same list as in Diraceq. 

open : opens a symbol, e.g. substitutes a combination of coupling constants like vinaen=ven-aen. 

AVALUE = (a) 

openall : opens all coupling constants, charges, etc. and then substitute s them, 
opensymbol : opens all symbols as openall but only in terms containing symbol 'a'. 
AVALUE = (a) 

p2D : replaces pi and pj by vectors Q and D: Q = —{pi + pj) and D = pi — pj. 

AVALUE = (i,j) 
p2l : changes a p^ to an invariant /. 

AVALUE = (p,I) 
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p2m : puts a 4-momciitum p to its mass shell, = —mp^, in the expressions and in the arguments of 
all functions. 

AVALUE = (p,mp) 

p2p : changes a p"^ to and p — > P in the string of gamma matrices. 

IVALUE, AVALUE = (I,p,P). 

1=0 p"^ changes to 

1=1 p2 changes to P^ and p ^ P 

p2Qs : expresses all scalar products Pi-pj in terms of pi, p^ and Q'^ = {p\ +P2f' for a three point function 
with Q +pi +P2 = 0. 

AVALUE = (p,mp) 

PoleSep : separates the PV functions explicitly into their 1/e pole parts and finite parts ^40^ and Bq . 
Pulling : is applied to expressions of the form 

with the result that p\ is placed next to u{p{) and P2 is placed in front of u{p2), after which the 
expression is simplified using the Dirac equation by a call to procedure Diraceq. 

IVALUE = (I) 

1=0 main option, is used in all programs up to ffbb boxes; 

eliminates j<4 in ii current and p2 in jj current 
1,2,3,4 is used in ffbb boxes; 

eliminates p\ or p2 or p^ or p4 in ii current. ^ 

PullingOrder : this procedure orders 7 strings containing p, 7^ and 7^ into ^^ifryv with one of the three 
factors possibly missing. Thus the surviving expressions are "f^p'-fv, 'y^p, 'y^jv, and pjv 

AVALUE = {p,i^,iy) 

Reduction : it has options 1=0, 1: if 1=0, then the user can perform a prereduction^ ^^by hand^; 
if 1=1, the standard prereduction is done automatically. 

After a prereduction is done, the reduction is performed on integrals of the form 

J dodid2dz 

where {1, g^, q^qv, ■ ■ .} means one of the expressions: scalar, vector, tensor ^ and di are given by 

do = q^+m\, 
di = {q+Pif+ml, 
d2 = {q + Pi + P2f + ml, 
d3_ = {q+Pi+P2+ Pif + ml. 

^See examples of assignment of current labels li in Figs. 13 and 14. 

^The prereduction consists of simplifications, such as the replacement of /{q^ + m?) by 1 — m? /{q^ + m?). 
^In the present version of SANC we have tensors of up to the 4th rank and N-point functions for N up to 4. 
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The N-point function, i.e. a one-loop diagram with N external legs, is defined by the following 
diagram: 

Pi i 




As a result the integrals (60) are replaced by linear combinations of Passarino-Veltman functions 
{PV)ij e {Aij, Bij, dj, Dij} where i = 1, 2, 3 or 4 for vector, 2nd, 3rd or 4th rank tensor, 
respectively, and j is a sequential number. 

IVALUE = (I) 

1=0 works without internal Prereduction 

1=1 the standard internal Prereduction is performed 

substitute : substitutes an argument "a" (charge, isospin or coupling constant of a particle). For example, 
to substitute the isospin of the particle 'typeID'=ll we call substitute(i3('typelD')) with the result 
i3('typelD')=l/2. 

AVALUE = (a) 

Scalarizing : expresses all (PV)ij functions in terms of scalars PVq. Option 1=0 acts differently for boxes 
and N=2,3 point functions; for 4-point functions the masses of the particles with momenta piand 
P2 are set equal to zero, for 2,3-point functions Scalarizing is exact in masses. Four digit options 
are applied only for 4-point functions as explained below. In all cases Scalarizing is exact in masses. 
Various options have been introduced to save CPU time. 

IVALUE = (I) 

1=1234 all masses are different 
1=1134 masses mi and m2 are equal 
1=1133 masses mi,m2 and ms,m4 arc equal 
1=1232 masses m2 and 1114 are equal (for boxes T5, T6) 
1=0 masses mi and m2 are equal to zero (for boxes) 
The latter option must be used also for self energies and vertices. 

Scalarizingdp : Scalarizing of dpy PV series, see Ref. [1] for definitions of d =dp functions. 
IVALUE = (I) 

1=0 scalarizing of dp series setting masses mi and m2 equal to zero 
1=1 scalarizing of dp series exact in all masses 
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ScalarizingProj : acts in ffbb boxes of topologies Tl, T3, T5 and T6. At first it projects the Input 
expression into a number of terms according to different structures with corresponding coefficients. 
Next it starts to apply the procedure Scalarizing to each term. The choice of index of Scalarizing(l) 
depends of the type of box topology. Finally, it forms Output expression by summing up all terms. 

AVALUE = (k3min,k4min,Topology,NameInput,p,mu,nu,pl,p2,I,NameOutput) 

k3min,k4min the indices, defining a cluster, see items Topologies T1,T3 
and Topologies T5,T6 

Topology number of the topology (1,3,5,6) 

Namelnput name of input expression 

p,mu,nu the same arguments as in PullingOrder 

pl,p2 fermionic 4-momenta defining basis of structures 

I index of procedure Scalarizing 

NameOutput name of output expression 

Sea I prod : calculates scalar products for 4-point function with 4-momenta satisfying pi +p2 +P3 +P4 = 0. 

AVALUE = (p) 

p scalar products pi ■ pj 

K scalar products ■ Kj, etc. 

Sing : In procedure Sing, the dimension n is set equal to 4 — £*, and then the PV functions, multiplied by 
£*, are analyzed: if a PV function has a pole, then the product e*xPV is replaced by its residue, 
and finally s* is set equal to zero. 

Symmetrize : symmetrizes Passarino-Veltman functions. Thus Bq is symmetrized using the symmetry 
property Bo(Q^, mi, 7712) = Bq{Q'^, 7712, mi). 

IVALUE = (I) 

1=0 symmetrizes Bq functions 

1=1 symmetrizes Bq and Co functions. 

Xil : sets all gauge parameters ^,^a,'Cz equal to one. These parameters are present in all intermediate 
contributions in gauge but cancel in gauge invariant physical observables. 

4.3 Feynman rules 

The SANC collection of Feynman rules is based on the Standard Model Lagrangian in gauge with three 

gauge fixing parameters ^, ^^i, and [1]. 

Propagators 

Every propagator should be multiplied by a factor of l/(27r)^ i . 

The propagator of a fermion / is a non-commuting function. It is defined by the following SANC command: 

pr(k,p,ii), ► 

/ 

where k is the field index, p is the fermionic 4-momentum and ii is the fermionic current label. 
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The vector boson propagator is a conirauting function: 



pr(k,/x,i/,p), jj, 1/ 

where k is the field index, fj,, v are the corresponding Lorentz indices and p is the bosonic 4-momentum. 
Vertices 

In the presently available class of diagrams, vertices are of three kinds: boson- fermion-fermion (bff), 
thrcc-boson (bbb) vertices and four-boson (bbbb) vertices. A vertex is a non-commuting function. Every 
vertex should be multiplied by a factor of (27r)'' i . 

bff vertices 

The SAMC command for this type of vertex is: 




vert(i,l,-j,a,ii), 

J 

where i, j and I are field indices, a is a Lorentz label and ii is a fermionic current index. The first field 
index refers to a boson; the other field indices refer to the incoming I and outgoing —j fermion fields. 
bbb vertices 

The SANC command for trilinear vector boson vertices is of the following form: 
vert(i,-j,l,a,/x,z/,(5,pi,p2) 




where i, j and I are boson field indices, a,^i,v are Lorentz labels and Q,pi,p2 are incoming momenta 
such that Q + Pi + P2 = 0. Significant is that the triplets of arguments, {i,-j,l}, {a, /x, i/} and {Q, pi, 
P2} are written in the same order according to the rule: ''from ingoing neutral to ingoing negative charge 
flow " , where the positive charge flow is shown by the arrows in the diagram. 

In the vertex diagrams involving a Higgs boson or scalar unphysical fields (Higgs-Kibble ghosts), the 
Lorentz indices /i are shown in brackets since they are dummy or silent indices, kept in the vert command 
for formal reasons, whereas these vertices do not depend on /x. 

In the book [1] all tri-linear bosonic vertices together with their Feynman rules, involving Higgs bosons, 
scalar unphysical fields (fp and <f>^ and Faddeev-Popov ghosts arc shown. Many of those diagrams carry 
silent (or dummy) Lorentz indices as discussed above in connection with bff vertices, and also dummy 
4-momenta. 

bbbb vertices 

The SANC command to define this class of vertices is (see the generic diagram) 
vert(i,j,k,l,a,/3,/i, u), 

j,l3 ^ k,n 
where i, j, fc, I arc field indices and a, /3, /i, v arc corresponding Lorentz indices. 

In the book [1] all quadri-linear bosonic vertices together with their Feynman rules are presented. 
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5 Processes, available in SANC v. 1.00 



In this section we briefly discuss the available in QED and EW branches processes. In this paper we do 
not discuss processes of QCD branch which are scarce. 

The Fig. 18 shows the fully open menu for "Processes" in the QED branch of SANC whose structure 
we briefly describe. The Figures 19 show all available 3-leg and 4-leg EW processes. We will not describe 
them in this part of the description since after explaining the QED branch, they may be easily interpreted. 
Moreover, the QED branch has mostly a pedagogical purpose. This is why it is worth devoting some 
time to it already in this first part mostly dealing with Precomputation. However, we will not describe 
here the structure of corresponding modules, leaving this for a second part of the SANC description. QED 
processes are presented by three classes: 1) a heavy photon decay; 2) e+e^ annihilation into a lepton 
pair (including Bhabha scattering); 3) Compton-like processes, i.e. e+e^ 77 or some other cross 
channel. Note that by our convention QED contains massless photons and three generations of leptons. 
We consider, nevertheless, the decay of a heavy photons for pedagogical reasons. 

When we arrive via a menu sequence, e.g. QED Processes ^ A ^ 11 Decay, we normally see 
three modules FF, HA and BR. Modules FF compute the scalar form factors of a given process. As was 
already stressed, they are channel independent, modulo a crossing transformation. Then a FORTRAN 
code to compute them can be automatically generated as described in Section 6. Modules HA compute 
channel dependent HAs. The channel is evident for examples 1) and 2); for example 3) we have for the 
time being HAs for the annihilation channel e~^e^ 77 and for the inverse channel 77 — > e+e^. 

Both in QED and EW branches, modules BR Root 
compute analytically the contributions due to ac- ^ [U^ Sanc 
companying Bremsstrahlung. In this connection, 
present implementation of Bremsstrahlung into 
SANC is not homogeneous. As a rule, for all 1 — > 2 
decays we have both Soft and Hard photon contri- 
butions. If only neutral particles are involved, for 
example Z — > i/v, the module BR is not present in 
the tree. There is one exception: for Z — > W^W~ 
decay we have only one module FF, since it is un- 
physical and we implemented FF for future use as 
a building block for more complicated processes. 
As a rule, for 4/ processes we also have both Soft 
and Hard photon Bremsstrahlung with exception 
of Bhabha process where we have only Soft con- 
tribution. For CC 2/1 2/ processes where we 
have realized quite involved calculations of Hard 
Bremsstrahlung with a possibility to impose sim- 
ple cuts, see also [29]. For the tree body decays 
t bli/ we have implemented both Soft and Hard 
contributions. For 2/2& processes we had so far no 
BR and no s2n calculations. 

This work is being started in version 1 . 10, 
where we have implemented FF and HA modules 
for three more ffbb processes: /i/i ZZ, /i/i Figure 18: Available processes in QED part. 

ZH and H j\j\Z . We recall, that /i stands for a massless fermion (its mass is retained only in the 
arguments of logs, if necessary). For the decay channel H — *■ fifiZ and annihilation process /i/i — *■ ZH 
we have implemented (BR) modules and s2n calculations, see also [30] which contains an extensive 
presentation of numerical results. 
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Root 
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Figure 19: Available processes in EW part. 
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6 User Guide 



6.1 Getting started 

6.1.1 SANC installation 

To work with SANC, one must install a SANC client on ones computer. The SANC client can be downloaded 
from the SANC project homepage http://sanc.jinr.ru or http://pcphsanc.cern.ch. On the homepage select 
Download, then download the client appropriate for your operational system (Linux, Windows), save it 
to your home directory and follow the instructions. ^ 

In Linux, opening the *.tgz file creates the directory /home/<user>/sancJnstaller. Change to that 
directory, read the Readme and execute install . sh. SANC is then installed in /home/<Mser>/sanc_1.00.00. 
To start a SANC session, go to directory /home/<Mser>/sanc_1.00.00/bin and execute sane. 

In Windows, start sanc_installer.exe program and follow the instructions, restart computer. To 
start a SANC session, click SANC client icon. 

6.1.2 SANC windows 

At the beginning of a client session the main SANC window opens, see Fig. 20,^ with several toolbars 
and windows or fields: 

• on top is the Menu bar with menus File, Edit, Build, Applications and View; 

• underneath is a row of three Toolbars: File, Edit and Build 

• underneath that on the left is the SANC tree field, and to the right of it the Editors List window; 

• underneath is the Output window and underneath that is the Console; 

• below, at the bottom, lies the Status bar. 
Other fields do arise in the course of the work. 

The five menus have the options shown in Table 1. Menus with — > have further extensions. For 
example. Toolbars has four options; they duplicate the File, Edit and Build toolbars, which are activated 
by default, and a latent option Memory. When the latter option is activated, two numbers are displayed: 
the first one is the current usage of the Java Virtual Machine (JVM) memory, and the second one is the 
total size of the JVM memory. All options can be unchecked in menu View Toolbars 

6.1.3 Login procedure 

• To log in, click the Login icon (the first icon of the File toolbar). The Login panel opens with a choice 
of SANC servers: local, sanc.jinr.ru and pcphsanc.cern.ch; choose one of the latter ones (the local server is 
for PCs which have the server itself installed), then enter the login name guest and password guest. 

• Click the Open Project icon (the second icon of the File toolbar). This opens the Open Project 
panel. There are two projects: Lessons and SANC. Select project SANC and press OK; then the SANC tree 
appears in the SANC tree field. 

6.1.4 The SANC tree 

The SANC tree has three options: QED, EW and QCD. Selection of one of these opens the next level 
of options: Precomputation and Processes. 

^To install and run SANC client one should have the Java Runtime Environment (JRE) at least version 5.0 Update 5 
installed, sec section Minimum System Requirements of the Download page at the SANC project homepage. 
®In the figure the windows are shown after several of the steps described below. 
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Tablo 1: Tlio SANC Menus and tlicnr options. 
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Here we describe the sequence of steps for option EW > Processes. The use of the Precomputation 
branch was described to an extent in Section 3. 

The available processes are subdivided into 31egs and 41egs. The two branches of 31egs are 3b and 
b2f decays, and those of 41egs are 4f and 2f2b processes; here b and f denote any boson and fermion, 
respectively. For each of the latter two there is a branch for Neutral Current and a branch for Charged 
Current processes. The next branching is into the available processes of that class. 

6.1.5 Naming conventions 

In SANC we use naming conventions for fields (or particles) shown in Table 2 where N is the field index, 
and in the columns headed "name" we show the names used internally in SANC. All associated parameter 
symbols are derived from these names. Thus the mass, charge and weak isospin of the electron are 
denoted mel, qel and i3el, respectively, also the vector and axial vector coupling constants (vel, ael) 
and their sum (vpael) and difference (vmael). 



Table 2: List of fields. 
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6.2 Benchmark case 1: b ff decays 
6.2.1 Semianalytical calculation 

Consider the Z ^ bb decay. First we open the relevant branch of the SANC tree: 

EW Processes b2f -> Z -> ff 

There are three FORM programs: (FF) Form Factor, (HA) Helicity Amplitudes, and (BR) Bremsstrah- 
lung. 

Select (FF) by a click with the right mouse button, this also pulls down a menu. On the menu left- 
click on Open. A Source Editor window opens with three tags: Form Editor, Fortran Editor, and 
Monte Ccirlo Editor. The first of these is activated by default and the FORM source code is displayed 
in the field. 

The particle indices can be seen in the Console field; by default they are: typeB = 2 {Z boson), 
typeU = 22 and typeD = 22 (& quarks). To change the final state fermions, their particle numbers can 
be changed by editing them in the Console field and pressing Enter. ^ 

After choosing the process, open the Numeric Form panel from the Application menu. In this panel 
the particle masses and other relevant information are displayed. 

Next the FORM code is compiled by clicking on the Compile button — the first icon in the Build 
toolbar (or by pressing the F7 function key). After compilation the FORM log file is shown in the 
Output field. 

Clicking on the Run S2N button generates the FORTRAN code; the FORTRAN code can be seen 
in the Output field. 

Repeat the sequence of steps for (HA) and (BR). 

The progress of work can be monitored by activating the Processes Table (see Table 1). 

The entire Output field is arranged in sheets with tags; for inspection any sheet can be brought to 
the forground by clicking on its tag. 

Once the three FORM codes (FF), (HA) and (BR) have been compiled and tranfered to the FOR- 
TRAN codes one can get the numeric results by the following sequence of operations: 

(i) open the FORTRAN editor sheet of the Editors List (belonging to the (FF) FORM code), 

(ii) open the Numeric Form panel from the Applications menu, 

(iii) press the Rehash button at the bottom of the Numeric Form panel, then the Compile button. 

The answer appears in the Output field. It consists of a list of Input parameters and a set of results: 
r(born), the total width [TotalWidth] in Born approximation, r(born-|-virt-|-soft) and the total width, 
r(born-|-virt+soft-|-hard). Also shown is the parameter w, set to 10~^°GeV by default. This parameter 
defines the separation between soft and hard radiation. It can be modified in the corresponding box of the 
Numeric Form panel. Rerunning the program after changing the value of ui (using the sequence Rehash 
> Compile) gives a result that differs only in the value of r(born+virt-|-soft). The born+virt+soft width 
is sensitive to parameter cu and can become unphysical (negative) for very small values of u). Increasing 
oj and rerunning gives positive values. 

^This need be done only once for a particular choice; to repeat, put the cursor at the end of the Console field, press the 
down-arrow key, select the required line of particle numbers using the up- and down-arrow keys and confirm by pressing 
Enter. 
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Ta1)l(- 3: Bouclmiark R(-sultK for TfZ ^ hb) dcaiy 





^ Born 


^ Born-\-virt-]-soft 


^Tolal 


SA 

MC 100 k 


0.356948 


0.336732 


0.360224 
0.360229 ± 0.000721 



6.2.2 Monte Carlo calculation 

The user can also carry out a Monte Carlo calculation generating various histograms: Photon Energy, 
Fermion Energy, Photon-Fermion Angular and Fermion-antiFermion Angular. To do this 

one must bring the Monte Carlo sheet of the Numeric Form to the foreground, check the boxes of 
the required histograms, then bring the MC sheet of the Editors List to the forground and rerun the 
program by choking the Compile button. After a while the Histogram Form is displayed. This form 
has a menu bar; menu Option allows display of the histogram statistics. On the Monte Carlo sheet one; 
can also select the random number generator,^ modify the number of MC events and the range of real 
photon energies fcomin and komax, where komin = and ho max can be used as an experimental cut. 
The Rehash button must be pressed after each change in the Numeric Form before clicking on the 
Compile button. 

The results for the decay rates (in GeV) of the semianalytical calculation and of the Monte Carlo 
calculation for 100000 events are summarised in Table 3. The numerical values are truncated to 6 
significant figures. 

6.3 Benchmark case 2: the process 2/ — > 2/ 

Consider the 4/ CC process — > //'. Implemented are the processes u d —> ^+ u^, its charge conjugate 

and the decay t ^ b£^ vg. For each process there arc three FORM programs: (FF) Form Factor, (HA) 
Helicity Amplitudes, and (BR) Bremsstrahlung. Each of these in turn is opened, compiled and run as 
above in Section 6.2. 

For process u J — > e+ i/g we have in the Console window the particle indices shown in Table 4. These 



Table 4: Assignment of particle numbers for process ud^ Vf. 



typelU = 


14 


initial Up-type antiparticle [d quark) 


typelD = 


13 


initial Down-type particle {u quark) 


typeFU = 


12 


final Up-type antiparticle (positron) 


typeFD = 


11 


final Down-type particle (neutrino) 



can be changed to typelU = 13, typelD = 14, typeFU = 11 and typeFD = 12 for process wrf — > e~ S-e by 
editing the particle numbers as explained above". 

Next bring the Fortran Editor sheet of the Editors List and the Numeric Form panel to the 
foreground. Shown on the Numeric Parameter sheet are the particle masses in GeV /c? and the CMS 
energy in GeV, also the cosine of the CMS angle between the incident and outgoing particle momenta. 

Click on the Rehash button at the bottom of the Numeric Form panel: the main module of 
FORTRAN code appears in the Fortran Editor sheet of the Editors List. Then click on Compile. The 

^Three random number generators provided arc; Ranlux, Ranmar and Mersenne Twister. 
^See Fig. 13a) for definitions of particle types typeXX. 
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final answer appears in the Output field. It consists of the parameters used (a, Gp, particle masses, the 
't Hooft scale /i and the Mandelstam variables), and the resulting differential cross sections da/d cos 
in picobarns in the Born approximation and Born+one-loop. The results for the default parameters and 
for several scattering angles are summarised in Table 5. The numerical values are truncated to 6 figures. 



Table 5: CMS differential cross sections in pb for ud — > e+fe 



cos 6 




ViGeV 
40 80 120 


-0.9 


Born 
Born -|- one-loop 


3.58202 10818.4 12.0561 
3.53427 9990.97 28.1226 


-0.5 


Born 
Born + one-loop 


2.23256 6742.78 7.51423 
2.18961 6226.00 12.8563 


0.0 


Born 

Born + one-loop 


0.99225 2996.79 3.33966 

0.97192 2769.11 5.12160 


0.5 


Born 
Born + one-loop 


0.24806 749.198 0.83491 
0.24453 695.224 1.47999 


0.9 


Born 
Born + one-loop 


0.00992 29.9679 0.03339 
0.01072 30.2878 0.09277 



Here the one-loop corrections are purely weak and QED corrections comprise one-loop virtual QED 
corrections and soft and hard radiations. 

The Born-|-one-loop cross section is insensitive to the 't Hooft scale parameter fj, which cancels between 
one-loop electroweak and the QED part of virtual corrections. 

Input parameters can be changed by editing the appropriate field of the Numeric Form panel and 
pressing the Rehash button. Again the Rehash button must be pressed before pressing Compile. 

In the NC sector there are many more processes. Here /i is a massless fermion of the first generation}^ 
or any neutrino, and / is any fermion. All procedures described above for the CC processes apply also 
in this case. 

Monte Carlo calculations are not yet implemented for 2 ^ 2 processes. 
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To produce whole Table 5 one can set flag tbprlnt = 1 in the Fortran Editors sheet. After editing the code one has 
not need press the Rehash button, but just Compile. 

'^^The masses of first generation fermions are retained only in logs to regulate coUineaj- singularities. 
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